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Abstract 

Neutrinos with energy greater than GeV are copiously produced in the p(A,p) interactions 
occurring in several astrophysical sites such as (i) the earth atmosphere, (ii) our galactic plane 
as well as in (iii) the galaxy clusters. A comparison of the tau and mu neutrino flux in the 
presence of neutrino oscillations from these three representative astrophysical sites is presented. It 
is pointed out that the non-atmospheric tau neutrino flux starts dominating over the downward 
going atmospheric tau neutrino flux for neutrino energy E as low as ~ 10 GeV. This energy value 
is much lower than the energy value, E > 5 x 10 4 GeV, estimated for the dominance of the non- 
atmospheric mu neutrino flux, in the presence of neutrino oscillations. Future prospects for possible 
observations of non-atmospheric tau neutrino flux are briefly mentioned. 

PACS numbers: 98.38-j, 13.85.Tp, 14.60.Pq 



* E-mail: athar@phys.cts.nthu.edu.tw 
t E-mail: cskim@yonsei.ac.kr 



1 



I. INTRODUCTION 



A present day main motivation for the extra-terrestrial neutrino astronomy is to obtain 
the first evidence of tau neutrinos from the cosmos around us above the relatively well known 
atmospheric neutrino background [lj. The tau neutrinos are essentially an unavoidable 
consequence of almost maximal neutrino flavor mixing between and v T as suggested by 
the atmospheric mu neutrino data analysis of the high statistics Super-Kamiokande detector 
(SKK) Q. 

A recent SKK analysis of the L/E distribution of the atmospheric mu neutrino data 
imply the following range of neutrino mixing parameters 

1.9 x 10~ 3 eV 2 < Am 2 < 3.0 x KT 3 eV 2 , sin 2 29 > 0.9. (1) 

This is a 90% CL. range with the best fit values approximately given by Am 2 = 2.4 x 
10~ 3 eV 2 and sin 29 = 1 respectively. This range of neutrino mixing parameters results 
in purely two flavor oscillation explanation of the L/E distribution of the atmospheric mu 
neutrino flux disappearance. The tau neutrinos as a result of these — > Vj- oscillations 
are so far identified on statistical basis (rather than on event by event basis) |4|. The total 
number of observed atmospheric non-tau neutrinos, on the other hand, are by now greater 
than 10 4 from various detectors ranging in energy between ~ 10 _1 GeV to ~ 10 3 GeV ^J. 
Given the recent detector developments, it is of some interest to estimate the tau neutrino 
flux from the earth atmosphere as well as from the nearby astrophysical sites in order to 
provide a more complete basis for the hypothesis of — > v T oscillations. 

The expectations for high energy astrophysical mu neutrino flux with E > 10 3 GeV with- 
out neutrino oscillation effects are summarized in 0. A purpose of this paper is to point 
out through examples that the GeV to TeV (1 TeV = 10 3 GeV) tau neutrino astronomy 
can be quite different as compared to the below TeV mu neutrino astronomy which is 
essentially dominated by the study of atmospheric neutrinos only. The recently discovered 
neutrino flavor oscillations does not bring any significant changes as far as the below TeV 
extra-terrestrial mu neutrino search/astronomy is concerned. Considering a different neu- 
trino flavor in the above energy range may change the situation somewhat by opening a 
new window to study cosmos. In this context, we shall distinguish between the possibilities 
offered by the different neutrino flavor astronomy through some examples. For a summary of 
above TeV astrophysical tau neutrino fluxes, see ^|. In this letter, we shall focus on the tau 



neutrinos and estimate the total tau neutrino flux from the three different astrophysical sites 

in order to illustrate that the GeV to TeV tau neutrino astronomy may be quite different 

from the GeV to TeV mu neutrino astronomy. 

The neutrino oscillation probability in the two neutrino flavor approximation is 

2 „ 2 / Am 2 (eV 2 )L(km)\ 
P(u, -> u T ) = sin 2 29 sin 2 (l27 ^ ^ j j . (2) 

Here L, is the neutrino flight length. In the earth atmosphere, it can be estimated using 



L = ^(h 2 + 2R®h) + {R® cos £) 2 - R® cos £. (3) 

The L is essentially the distance between the detector and the height at which the atmo- 
spheric mu neutrinos are produced. The Rq ~ 6.4 • 10 3 km is the earth radius, and h — 15 
km is the mean altitude at which the atmospheric mu neutrinos are produced. In general, h 
is not only a function of the zenith angle £, the neutrino flavor but also the neutrino energy 

a. 

This letter is organized as follows. In section 2, the mu and tau neutrino flux originating 
from the earth atmosphere, the galactic plane as well as from the galaxy clusters is briefly 
discussed. In section 3, the neutrino oscillation effects are studied for these. In section 4, 
the limited future prospects for possible observations of non-atmospheric tau neutrinos are 
mentioned, whereas in section 5, conclusions are presented. 



II. SOME EXAMPLES OF GEV TO TEV ASTROPHYSICAL TAU NEUTRINOS 

In order to estimate the total tau neutrino flux from a given astrophysical site which is 
composed of intrinsic and oscillated ones, we need to first estimate the intrinsic mu and tau 
neutrino flux from that site. In this section, we briefly describe the general procedure used 
to estimate these. 

Before doing that, let us remark here that the incident cosmic ray energy that we consider 
here ranges typically between 10 < E p /GeV < 10 4 . The corresponding cm energy range 
in pp collisions is 1 < y / s^/GeV < 10 2 as s pp ~ 2m p E p . The data from various collider 
experiments indicates that the tau neutrino parent hadron (Dg) production cross section 
here is suppressed relative to mu neutrino parent hadron (tt^) production cross section in 
pp collisions in the above yfs^ range [lOj. Therefore, any sizeable tau neutrino flux from 
the following sites should be verifying the neutrino flavor mixing. 
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A. Atmospheric 



Briefly, the incoming cosmic rays interact with the air nuclei A, in the earth atmosphere 
and give rise to mu neutrino flux. For 1 < E/GeV < 10 3 , the n^, K production and their 
direct and indirect decays are the main sources of mu neutrinos, both being in the region of 
conventional mu neutrino production 11]. The absolute normalization of the conventional 
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atmospheric neutrino flux is presently known to be no better than (20—25)% 

For present estimates, the mu neutrino flux is taken from Ref. jl3(- These are neutrino 
flux calculations in one dimension without geomagnetic field effects. The up down mu 
neutrino flux is taken to be the same, as the present discussion is independent of any specific 
detector. At higher energy, the prompt mu neutrino production from D's dominates over 
the conventional one \\\ - 

The atmospheric tau neutrino flux arises mainly from and is calculated in Ref. 15. 
The Quark Gluon String Model (QGSM) is used in Ref. [l6| to model the pA interactions. 
The low energy atmospheric tau neutrino flux is essentially isotropic 15]. For E < 10 3 
GeV, the atmospheric tau neutrino flux is obtained by following the procedure given in 



Ref. 
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161 ] and re-scaling w.r.t new cosmic ray flux spectrum, taking it to be dominantly 



the protons jl8, 
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20]. 



B. Galactic Plane 



The galactic plane mu neutrino flux is calculated in Ref. [21|, |22j , whereas the galactic 
plane tau neutrino flux is calculated in Ref. Q] . These calculations consider pp interactions 
inside the galaxy with target proton number density ~ 1/cm 3 along the galactic plane, 
under the assumption that the cosmic ray flux spectrum in the Galaxy is constant at its 
locally observed value. The current Energetic Gamma Ray Experiment Telescope (EGRET) 



observations of the diffuse gamma-ray emission from the galactic p 



ane seems to imply a 



3. 



We have checked 



slightly less (~ 0.62) target proton density along the galactic plane 
that taking into account this in our estimates makes only a minor difference. 

Following Ref. Q], the galactic plane mu and tau neutrino flux for E < 10 3 GeV is 
obtained by re-scaling w.r.t new cosmic ray flux spectrum. The tau neutrino production is 
rather suppressed in the galactic plane relative to mu neutrino production. The mu neutrino 
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flux is larger than the tau neutrino flux for E < 10 3 GeV from the above two sites. 

Figffl gives the intrinsic mu and tau neutrino flux, F V (E) = diV°/d(logio-E) in units of 
cm~ 2 s _1 sr _1 , estimated using the above description. The figure also shows the cosmic-ray 
proton flux spectrum we have used. 



C. Clusters of galaxies 

Clusters of galaxies are presently considered to be an interesting laboratory to investigate 
various sta ges of structure formation via the study of diffuse gamma-rays (and cosmic-rays) 
from these j^H]. 

The mu and tau neutrino fluxes are produced here in pp interactions (of cosmic-rays with 
the intra cluster gas), under the main assumption that a large fraction of the cosmologically 
produced baryons is inside these galaxy clusters. The intrinsic tau neutrino flux estimate 
here is rather similar to the galactic plane situation. 

The maximum mu neutrino flux from clusters of galaxies is estimated by correlating it 
with the corresponding gamma-ray flux in pp collisions. This estimate is under the assump- 
tion that diffuse gamma-ray flux observed in the energy range 10 < E/GeV < 50 by EGRET 
is dominantly from clusters of galaxies. We shall use this upper limit as an another example 
of an astrophysical site that may produce tau neutrino flux in the GeV to TeV energy range. 

The detailed model dependent calculations indicate that the actual contribution of the 
galaxy clusters towards the recently observed extra-galactic diffuse gamma-ray flux is ap- 
proximately two orders of magnitude small j^ . 

Before studying the effects of neutrino oscillations on these intrinsic neutrino fluxes, let us 
here emphasize that the non-atmospheric mu and tau neutrino fluxes can only be considered 
as upper limits on these fluxes in light of the existing inherent uncertainties in estimating 
these. 



III. EFFECTS OF NEUTRINO OSCILLATIONS 



We shall perform here the two neutrino flavor oscillation analysis. In the context of two 
neutrino flavor oscillations, there are only two neutrino mixing parameters. The mixing 
angle 9 and the Am 2 (= m\ — m\). There are no matter effects for these two flavors. 
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In the two flavor approximation, the total tau neutrino flux from an astrophysical site 
is given by 



where P^E) = P(z/ M -> v T ) is given by Eq. © and P(z/ M -»• u^) = 1 - P(z^ -> i/ r ). 
The F®(E) is the intrinsic neutrino flux in units of cm -2 s~ 1 sr~ 1 and is taken according to 
discussion in section 2. 

Three general directions in the earth atmosphere are considered as representative ex- 
amples to compare the atmospheric tau neutrino flux with the non-atmospheric one in the 
presence of neutrino oscillations. These are the downward, the horizontal and the upward 
directions. Fig. El depicts the L osc given by Eq. (J2J) for the range of Am 2 given by Eq. (JTJ) 
with maximal mixing. The three distances are taken from Eq. @5J), with, for instance, the 
downward distance is obtained by setting £ = 0. The horizontal distance is obtained by 
setting £ = 7r/2. 

Using Eq. (flj). the total downward going atmospheric tau neutrino flux is estimated. It 
is then compared with the total galactic plane and total upper limit galaxy clusters tau 
neutrino flux in Fig. El for the whole range of Am 2 with maximal mixing. The distance L 
for galactic plane neutrinos is taken as ~ 5 kpc, where 1 pc ~ 3 x 10 13 km. For galaxy 
clusters, its representative value is taken as 1 Mpc. Since L osc <C L, the galactic plane and 
galaxy cluster mu neutrinos oscillate before reaching the earth. Also, note that this flux is 
averaged out for the whole range of Am 2 in the entire considered energy range. The effect 
of different Am 2 values diminishes for E > 50 GeV for total atmospheric tau neutrino flux. 
From the figure, it can be seen that the galactic plane/non-atmospheric tau neutrino flux 
starts dominating over the downward going atmospheric tau neutrino flux even for E as 
low as 10 GeV in the presence of neutrino oscillations. This is a very specific feature of 
tau neutrinos, and is absent for mu neutrinos. This specific behavior has to do with the 
neutrino oscillations. The galactic plane tau neutrino flux for 1 < E/GeV < 10 3 in the 
presence of neutrino oscillations can be parameterized as 



where F Vr is in units of cm~ 2 s _1 sr _1 and on r.h.s. E is in units of GeV. 

In Fig. 01 the galactic plane tau neutrino flux is compared with the atmospheric tau 
neutrino flux, using Eq. ^ for the three general directions for the atmospheric tau neutrino 



F Vt {E) = P^E) ■ F«{E) + P m (E) ■ P° (P) 



(4) 



F Ut {E) = 1.31 • 10" 5 ■ E im P + 2.15exp(-0.21\/P) 



1 -2.74 



(5) 
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flux reaching the detector. Here the best fit values of the neutrino mixing parameters are 
used. The oscillatory nature of the upward going tau neutrino flux can be seen from Eq. 
(J2J). The cross over for the galactic tau neutrinos relative to the horizontal atmospheric 
tau neutrinos occurs at ~ 50 GeV, whereas the same occurs for the upward direction at 
~ 400 GeV. The total atmospheric tau neutrino flux is maximum in the upward direction. 
It is minimum in downward direction, relative to the galactic plane tau neutrino flux in the 
presence of neutrino oscillations, owing to the behavior of L/L osc ratio as a function of 
neutrino energy. Fig. 0] indicates that zenith angle dependence of the total tau neutrino 
flux can at least in principle help to distinguish between atmospheric and non-atmospheric 
tau neutrino flux. The galactic tau neutrino flux transverse to the galactic plane is three to 
four orders of magnitude smaller than the galactic plane one jlq . 

Fig. El gives a comparison of the downward going atmospheric and the galactic plane 
mu neutrino flux in the presence of neutrino oscillations. For this comparison, mu neutrino 
flux is taken from Ref. [3] without re-scaling for E < 10 3 GeV. This mu neutrino flux 
includes contribution from the .D's for E > 6.3 x 10 5 GeV. The total mu neutrino flux is 
estimated according to Eq. (jlj) with appropriate modifications for the best fit values of 
the two neutrino mixing parameters. In contrast to the possibility of seeing the galactic 
plane with multi GeV tau neutrinos, note here that with mu neutrinos, it can occur only 
for E > 10 5 GeV. 

A relevant remark is that for the best fit values of the neutrino mixing parameters, 
the Piy^ — > v T ) is relatively large along the horizontal and upward directions in the earth 
atmosphere [see Eq. (j2J) and Fig. ^ for 1 < E/GeV < 10. So essentially the atmospheric mu 
neutrino flux in the absence of neutrino oscillations alone determine the total atmospheric 
tau neutrino flux in comparison with the total galactic tau neutrino flux. 



IV. PROSPECTS FOR POSSIBLE FUTURE OBSERVATIONS 

In this section, we shall estimate the galactic plane tau neutrino induced shower event 
rate for a megaton class of detectors to indicate the limited prospects offered by GeV to 
TeV tau neutrino astronomy to search for extra atmospheric astrophysical neutrino sources 
in this energy range. 

Let us add a remark here that, at present, the dedicated large high energy neutrino 
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detectors such as the Antarctic Muon And Neutrino detector Array (AMANDA) are also 
sensitive to all neutrino flavors, however, with energy ~ 50 TeV [^\ - 

For 10 < E/GeV < 10 3 , a signature for the tau neutrinos is to measure the energy 
spectrum of the tau lepton induced electromagnetic and hadronic showers produced in tau 
neutrino nucleon interactions occurring inside a densely instrumented Cherenkov radiation 
detector 



j. Though, it is a challenging task to distinguish between tau and non-tau 
neutrinos for the present generation of detectors in the above energy range J^, however 
certain shower signatures remain distinctive for tau neutrinos 2s| |. 

The galactic tau neutrino induced shower production rate can be approximately estimated 
by convolving the total galactic tau neutrino flux in the presence of neutrino oscillations, 
given by Eq. © with the a^ N (E). The a^ N (E) for 10 < E/GeV < 10 3 is parameterized 
as 

a™ N {E) = -4.43 + 0.52 • E + 3.58 • 10" 4 • E 2 . (6) 

The <J® N is in units of 10 -38 cm 2 and E is in units of GeV. The CTEQ6 parton distribution 
functions 



30| were used to estimate the cross section. 
For recent detailed evaluations of cr„ c N , see Ref. j3l|. The possible tau lepton polarization 
effects (32! are not taken into account in the event rate estimates presented here. 

Table I gives the galactic tau neutrino induced shower event rate for a one Mega ton 
detector, in units of (Mt -yr) _1 in 27r steradians of upper hemisphere in eight logarithmically 
equally spaced energy bins. The table indicates that the per year event rate is about 1.5. 
With a 3 to 5 year data collection time for a one Mega ton detector, the galactic tau 
neutrino induced shower event rate can thus be in the range of ~ O(10) for E > 10 GeV. 
This detector faces only the downward going atmospheric tau neutrino flux as background 
to the dominant galactic plane tau neutrino flux in the presence of neutrino oscillations. 



V. CONCLUSIONS 

1. The effects of neutrino oscillations on low energy (1 GeV < E < 10 3 GeV) tau 
neutrino flux produced in the earth atmosphere, in our galactic plane and in galaxy clusters 
are presented in two neutrino flavor approximation. 

2. The galactic plane should be observable with tau neutrinos with energy > 10 GeV, 
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depending on the orientation of the concerned detector w.r.t. galactic center/plane at the 
time of observation. The observation of galactic plane with multi GeV tau neutrinos is in 
sharp contrast to the case of mu neutrinos with which the galactic plane is observable only 
with energy > 10 5 GeV for the same orientation of the detector. 

3. Transverse to the galactic plane, the maximum galaxy clusters tau neutrino flux is also 
above the atmospheric one, providing another example of a new window to study cosmos in 
the above energy range. 

4. This observation may also have some relevance for the long baseline experiments 
searching for the tau neutrinos in — ■> v T oscillations j^ . 
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TABLE I: Galactic plane tau neutrino induced shower event rate for a Megaton class of detectors 
in eight logarithmically equal energy bins. Details are given in the text. 



Energy Bin (GeV) 


Event Rate (Mt • yr • 27rsr)" 1 


Energy Bin (GeV) 


Event Rate (Mt • yr • 27rsr) _1 


10 - 17.78 


0.22 


100 - 177.83 


0.16 


17.78 - 31.62 


0.30 


177.83- 316.23 


0.12 


31.62 - 56.23 


0.26 


316.23- 562.34 


0.10 


56.23 - 100 


0.22 


562.34-1000 


0.07 
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Intrinsic Galactic Plane Neutrino Flux 
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FIG. 1: The intrinsic galactic plane mu and tau neutrino fluxes. The cosmic-ray flux spectrum 
used in the estimates is also shown. 



13 



E 



1e+08 
1e+07 r 
1e+06 r 



~ 1e+05 r 



1e+04 r 



1e+01 



-i 1 1 1 — i — i — i — r 



n rzz — i 1 — i — i — i — r 



L (Am =1.9-10 3 eV 2 ) 
Downward distance 
Horizontal distance 
Upward distance 



-cr 1 1 1 — i — i — i — r~ 




_i i i i i i i i_ 



_i i i i i i i i_ 



_i i i i i i i i_ 



10 



100 



1000 



E (GeV) 



FIG. 2: The -> v T oscillation length L osc [= 0.787 km £'(GeV)/Am 2 (eV 2 )] in km as a function 
of neutrino energy in GeV. The three general distances traversed by mu neutrinos in the earth 
atmosphere are also shown (the horizontal lines). More details are given in the text. 
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FIG. 3: An illustrative comparison of the galactic plane, the maximum galaxy clusters and the 
downward going atmospheric tau neutrino flux in the presence of neutrino oscillations as a function 
of neutrino energy. 
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FIG. 4: The galactic plane and the atmospheric tau neutrino flux in the presence of neutrino 
oscillations for the three general directions in the earth atmosphere as a function of neutrino 
energy. The neutrino mixing parameter values used here are the approximate best fit values, i.e., 
Am 2 = 2.4 x 1(T 3 eV 2 and sin 2 26 = 1. 
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FIG. 5: Comparison of the downward going atmospheric mu neutrino flux, the galactic plane 
mu neutrino flux and the maximum galaxy cluster mu neutrino flux in the presence of neutrino 
oscillations. The non-atmospheric mu neutrino flux starts dominating over the atmospheric one 
only for E > 5 x 10 4 GeV. 
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